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Abstract 

The spherical limit of strongly commensurate dirty bosons is studied perturbatively 
at weak disorder and numerically at strong disorder in two dimensions (2D). We ar- 
gue that disorder is not perfectly screened by interactions, and consequently that the 
ground state in the effective Anderson localisation problem always remains localised. 
As a result there is only a gapped Mott insulator phase in the theory. Comparisons 
with other studies and the parallel with disordered fermions in 2D are discussed. We 
conjecture that while for the physical cases N = 2 (XY) and N = 1 (Ising) the theory 
should have the ordered phase, it may not for N = 3 (Heisenberg). 



1 Introduction 

The problem of interacting bosons in a random potential is a paradigmatic case of an inter- 
acting disordered system, and as such has attracted much attention throughout the years 
0, 0. Although in one version or another it has been used to describe numerous physical 
situations || , it has proven very difficult for a theoretical analysis, since it inextricably com- 
bines the effects of interactions and Anderson localisation. Just like its fermionic cousin the 
metal-insulator transition the problem of dirty bosons seems to lack a simple analytic 
mean-field theory around which to begin a systematic study. Most of the information on the 
dirty boson quantum phase transitions derive therefore from numerical studies ||, and more 
recently from an expansion around the lower critical dimension |J. 

In this paper we will be concerned with a limited class of the dirty boson models at a 
commensurate filling, and study the limit where the number of bosonic species N is large 
0. As is well known in this limit the mean-field theory, or the saddle-point approximation, 
becomes the exact solution. The quantum mechanical action at T = that defines our 
problem is 

Sm = J d D xdr | (d T V(x , r)) 2 + (W(f , r)) 2 + (V(x) - /i) V 2 (x , r) + ^ 4 (x, r) j , (1) 
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N 

where *&(x, r) is a real A^-component bosonic field, \l/ 2 = ^ \I/ 2 , and V(x) is a random (in 

a=l 

space) external potential. For simplicity, it will be assumed that V(x) is uncorrelated, so 
that (V(x)V(y)) = W5(x — y). We will mostly be interested in two dimensions (D = 2), 
but will leave a general D in the action to comment later on results in other dimensions. 
Note that disorder is assumed to be a random function only of spatial coordinates, while it 
is completely correlated in (i.e. independent of) the imaginary time. This is what makes it 
much stronger than in the corresponding problem in classical mechanics. The theory (1) for 
N = 2 describes the superfluid order parameter in the Bose-Hubbard model, at a density 
of bosons commensurate with the lattice @, |§, also known in literature as the random-rod 
problem ||. For N = 3 the theory may be used to describe disordered quantum rotors, i.e. 
the magnetic quantum phase transitions in the Heisenberg universality class in the presence 
of quenched randomness ||10|| . When N — 1 the theory describes a random system with the 
Ising symmetry. In general the action (1) provides a minimal description of the quantum 
disordered interacting system, and for N = oo has been studied by renormalisation group 
methods in the past [JTTJ, [L2| with conflicting results. The purpose of this paper is to shed 
some light on the physics implicit in this model, and, in particular, to argue that the model 
allows no superfluid phase in D = 2. 

To see what is involved in solving the problem in the (spherical) limit N = oo, perform 
the standard Hubbard-Stratonovich transformation on the quartic term and integrate out 
all but one of the bosonic fields. This leaves one with the transformed action: 

S eS [ X ,i>] = J^xdr^-^i^^ + id^^r)) 2 (2) 
+(V*i(£, r)) 2 + (V(x) + X (x, r) r)j 
+ ^(iV - 1) lndet {-d 2 T - V 2 + V(x) + X (x, r) - //} , 

which is just the original problem rewritten exactly. Assuming that the Hubbard-Stratonovich 
field at the saddle-point is independent of imaginary time, x(x, r) = X (x), and that ^i(x, r) = 
N2C(f>o(x), the saddle-point equations become 



x{x) = A ( x , r 



f,r)+c0 o (f), (3) 



,Q2 _ v 2 + V(x) + X (x) - n 

e c = 0, (4) 

where 4> a (x) are the random eigenstates, and e a the random eigenvalues of the susceptibility 
matrix 

(-V 2 + V{x) + X(x) - /i) 4>a{x) = £ a <f> a (x), (5) 

with Eq being the lowest eigenvalue. The Eqs. ©-© are completely standard, and the 
only novelty compared to the case without disorder |13[ is the random spectrum instead of 
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the usual plane waves. In principle, one may expect to have two phases: sq ^ and c = 0, 
corresponding to the gapped Mott insulator (MI), or e = and c ^ which would represent 
a superfluid (SF). The gapless insulating Bose-glass (BG) between the MI and the SF 
should in general be absent, as we argue below. 

With V(x) = 0, the solution x(%) = Xo is uniform, and the model leads to the well- 



known large-N critical behavior in D + 1 dimensions [13[]. The correlation length exponent, 
for example, in the pure case is v = 1/(D — 1), and in D = 2 the Harris criterion |14j (that 
says that disorder is irrelevant if vD > 2) implies that disorder is precisely marginal. When 
V(x) 7^ 0, in the MI phase the saddle-point Eq. (||D, after integration over the frequency, 
can be written in the basis {</>«} as: 

XW = AL — ^=-- (6) 



The functions 4> a are the eigenstates of the screened, but nevertheless random, potential 
V{x) + x(x)-> an d would therefore naively all be expected to be localised in D = 2 [|15| . 
In particular, for the localised ground state the first term in the sum in Eq. (|6|) becomes 
large as 6q — > precisely in the region of localisation, which by self-consistency implies that 
x{%) is also large there. That, on the other hand, then implies Eq is large, and not small as 
assumed, and one runs into a contradiction. Evidently, for the spectrum to extend all the 
way to zero the discrete sum in the last equation must be approximable by an integral, so 
that the infrared singularity becomes integrable. For this to occur the weight of each of the 
terms corresponding to the low-energy states in Eq. (^) must vanish in the thermodynamic 
limit as the inverse of the system size, which is tantamount to delocalisation of the low- 
energy eigenstates. Put differently, the collapse of the gap must be accompanied by the 
simultaneous delocalisation of the ground state, so that the gapless phase is necessarily a 
SF. There can be no intermediate localised BG in the model at N = oo. 

With this picture in mind the appearance of the superfluid phase in the large-N model 
in D = 2 appears rather counterintuitive: although screening introduces correlations into 
the effective random potential, the states should nevertheless always remain localised. In 
the rest of the paper we first show that although to the lowest order screening does reduce 
the random potential, it does not make it completely smooth and consequently the MI gap 
can not close. This conclusion is further corroborated by the numerical solution of the self- 
consistent equations on a lattice and absence of the finite-size scaling of the gap and the 
ground state participation ratio. In the closing section we compare our result with other 
studies and speculate on the implications for physical cases N = 1,2, 3. 



2 Weak-disorder expansion 

For a given random configuration the self-consistent equations can not be solved analytically, 
and one has to resort to numerical computations. For weak disorder, however, we can 
expand the matrix element in (S) in powers of the screened potential. To that end write 
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x{%) = Xo + Xi(z)i where / xi(x)d D x = 0. The uniform part xo is just the renormalisation 
of the chemical potential, while V(x) = V(x) + Xi(af) is the screened potential, which should 
vanish with vanishing randomness. Expanding the right hand side of (H) in the MI phase in 
V(x) and taking the Fourier transform, we get (for q 0) 

V(q) = V(fi-\n(fiV(fi+\[dkh(k,q)V(k)V(q-k) (7) 
-A J dkdll 2 (k, I, q)V(k)V{l)V{q- k-l) + 0{V A ), 

where 

U(q) = J dpdcoG (u,p)G (u,p + q), (8) 
is the standard polarisation bubble, and 

Ii(k,q) = JdpduG (uj,p)G (uj,p + k)Go(uj,p + q), (9) 

and 

I 2 (k,l,q)= dpduGo(u,p^G (uj,p + k)Go(u,p + k + l)Go(u,p + q). (10) 

The propagator for the clean case is given by Gq 1 (uj,p) = uo 2 +p 2 + il 2 , where fi 2 = Xo—^ > 
and is the MI gap. Eq. (0) can be represented diagrammatically as in Fig. |l|. 

We next introduce the two point correlator W(q)5(r) = (V(q)V(—q + r}\ where (• • •) 
represents disorder averaging, as a measure of the screened disorder. From Eq. ([?]) it follows 
that 

W(q){l + XU(q)} 2 = W(q^+2x[dkI 1 {k,q^(V{-q)V(k)V(q-k)) (11) 

-2A J dMl 2 {kJ,q^(V{-q^V{k)V{t)V{q-k-V)) 
+A 2 J dkdk'h{k, q)h{k', -q) (V{k)V{q- k)V{k')V{-q- k' 



+0(W 3 ). 

Diagrammatically, the second-order contributions may be represented as in Fig. |2|. In 
the Appendix we compute the above averages in D = 2. Note that although the random 
potential is assumed uncorrelated in space, the screened potential develops correlations and 
W(q) becomes a non-trivial function of the wave-vector. For the low-energy states one 
expects the localisation properties to be determined by W(q) at small q, so we focus on the 
limit q — > and denote W(q — > 0) = W. To the second order in W in the limit Q — > and 
in D = 2 one then finds (see the Appendix for details): 

W ^2 ( W \ 2 n2 C 1 ^/A\ 2 32 /A 



4(£HK£))}+°<^ 
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where the constant c = 1/(8-71"), and A is the ultraviolet cutoff implicit in (0). 

The last equation is our central result, and several remarks are in order. To the first 
order in W, one finds that as Q — > 0, W — > 0, which one may be tempted to interpret as a sign 
of delocalisation of the ground state. This is a consequence of the screening of the random 
potential by the medium, which to the zeroth order in disorder is pure and thus screens 
perfectly at q = 0. Also, recognising the combination W/Q 2 as a dimensionless measure of 
screened disorder, to the lowest order Eq. (|i~2|) agrees with the Harris criterion: disorder is 
marginal in D = 2. The fate of disorder is therefore determined by the higher-order terms 
in the expansion. To the second order in disorder we find that 

~ A 2 / W \ 2 

W — » — =■ ( -To - ? ) , as Q^O, (13) 



2tt 5 U 2 c 2 

i.e. goes to a non-universal finite constant as the gap decreases. If the bare disorder is 
weak the screened disorder will be even weaker, but always finite. The consequence is that 
the ground state and the excited states in D = 2 should remain localised [jn|, so that our 



qualitative argument from the introduction would imply that the gap can not close. This is 
in accordance with the direct numerical solution at strong disorder to which we turn next. 



3 Numerical solution 

We begin by introducing a discrete version of our theory where the continuous variable x is 
replaced by a lattice-site index i on a quadratic lattice of linear size L. The kinetic energy 
term V 2 is replaced by the nearest-neighbour hopping measured by t, the random potential 
is chosen from a uniform distribution of width W and the interaction strength is given by A. 
In our calculations we set W/t = 4 and X/t = 8, which corresponds to strong disorder and 
interactions. After the integration over frequency, the self-consistent Eq. becomes 

,= A^, (14) 



a=l 

in the {(f> a (i)} basis where these wave functions are now eigenvectors of the matrix 

X! {'hi + ( V i +Xi~P) <f>a{j) = £a(f)a(i) (15) 

3 

where is non-zero for nearest-neighbour i,j only. 

We solve the set of L x L equations using the Newton- Raphson algorithm. We gradually 
increase the chemical potential //, using the last found solution as the initial guess at the 
next jji,. Finally, we average over many disorder configurations. Of course, for finite L the 
gap is always finite, so to infer the result in the thermodynamic limit we make the standard 
finite-size scaling ansatz for the average ground state energy 

e = L- z F\L$(fM-iM c j\, (16) 
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where z is the dynamical critical exponent, v is the correlation length exponent, and /z c is the 
critical point in the thermodynamic limit; F(x) is a universal scaling function. The values 
of z and /i c are determined by scaling the £o-axis until all curves cross at a single point. The 
exponent v is found by scaling the /i-axis so that a reasonable collapse of all the data onto 
a single curve is achieved. 

Such an attempt of finite size scaling of our data is shown in Fig. [3] for systems of linear 
size L = 6,8, 10, 12. We display the result for the value z = 0.9, but the picture remains 
qualitatively the same for all 0.5 < z < 1.0. The gap continuously decreases with /i, but the 
failure of the finite size scaling suggests it does not vanish in the thermodynamic limit. 

We have also argued that at the point of collapse of the gap, the ground state would 
be expected to become delocalised. A useful measure of the degree of localisation of the 
wavefunctions at a given energy is provided by the participation ratio 

P(s) = E 5 % a ~ £) , (17) 



H)\ 4 



i=i 



which is proportional to 1/L 2 for the localised states and approaches a constant for the 
extended ones. In the critical region, one expects the participation ratio to assume a similar 
finite-size scaling form: 



P(e ) = L-(°- D f^ L-^ii-ii, 
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where Df is the fractal dimension of the ground state wavefunction and <&(x) another scaling 
function. Our data for the participation ratio are shown in Fig. [| for the sizes L = 8, 10, 12. 
Again, attempts to find the common crossing point by tuning Df fail. We see that the 
participation ratio of the ground state grows as \i is increased, but conclude that the ground 
state nevertheless seems to remain localised. This is consistent with the data for the ground 
state energy. 

Our Newton-Raphson algorithm has difficulties converging as \i is increased and the 
problem becomes more non-linear. It is possible we simply have not been able to reach the 
critical point in our numerical calculation. When taken together with the weak-disorder 
expansion and the physical arguments, however, we believe a more likely interpretation is 
that there is no SF phase in the model. 



4 Conclusion 

To summarise, we studied the large-N limit of the commensurate dirty boson theory, and ar- 
gued that at weak disorder screening does not delocalise the ground state, and consequently, 
that there is no MI-SF transition in D = 2. Numerical results for the ground state energy 
and the participation ratio that support this conclusion were provided. 

Our conclusion agrees with the results of Kim and Wen [TJJ who found that disorder 



is always relevant for D > 2 and could not find any stable critical points within their 
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renormalisation scheme. The latter point may in principle be interpreted in three ways: as 
a failure of the renormalisation procedure, as that the transition is discontinuous, or that 
there is no transition. Our findings support the third conclusion. On the other hand, we are 
in disagreement with the recent study of Hastings who considered the closely related 



random spherical model, and found the disorder to be marginally irrelevant in D = 2. At 
the moment we do not fully understand what is the resolution of this disagreement, nor how 
the ground state becomes extended in Hastings' theory. 

While we were mostly concerned with D = 2, the same perturbative procedure can be 
repeated in D = 3. We found that the same diagram in Fig. |2|(e) that led to the finite term 
for W in D = 3 vanishes logarithmically as the gap decreases. More importantly, in D = 3 
the Anderson localisation problem allows a mobility edge, so the screened disorder need not 
go all the way to zero for the ground state to delocalise. We would therefore expect that the 
theory (1) would have a MI-SF transition in D = 3, as apparently has been found in earlier 
numerical calculations [M. 



An important question is what our considerations imply for the physical cases N = 1,2,3 
mentioned in the introduction. We believe that in D = 2, for N = 2 the theory (1) does 
have a transition and which is in the BG-SF universality class. This has been found in 
the dual theory for the commensurate dirty-bosons |J, in both D = 1 and D = 2, and in 
detailed numerical calculations || [L7|. The BG-SF transition is best understood in terms 
of disorder-induced proliferation of topological defects, and thus is very specific to having a 
complex (N = 2) order parameter. The same topological mechanism will not apply to the 
case of a random quantum ferromagnet N = 3, and we conjecture that for N = 3 there may 
not be a gapless phase in D = 2. On the same grounds, we expect that for the Ising case 
N — 1 the transition again will exist fTB |. 



Finally, we note the similarity between our problem and the problem of interacting 
disordered fermions in D = 2 |fT9| . In the large-N limit the metallic phase in the fermionic 
problem would correspond to an extended state at the Fermi level, as opposed to the extended 
ground state in our problem. Nevertheless, one can show [2(J that already to the lowest order 
in disorder, screened disorder remains finite, and thus the state should remain localised. We 
would therefore expect that the fermionic version of the action (1) also should have only the 
localised phase in D = 2, at least in the large-N limit. 
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6 Appendix 

In this appendix, we provide the details of the calculations leading up to our main analytic 
result (|I2"D. We begin by calculating the integrals (|§l)-(|T0l). Using the standard Feynman 
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parameters [[13] the integrals can be rewritten as 



(47r)^ [ t (! _ t ) q 2 + q 2 ] 

i["(f)-'(;) ! +°((ff 



dt- 



3-D 
2 



_c 



(19) 



A(M) 



i r 



2(4^ 



1 



c 
c 

4Q 3 ' 



dt 

[t(l - t)fc 2 + ft 2 ] 

4(f) 2 -16(f) 4 + 0((f 



5-g 
2 



k 



• oo 
0. 



(20) 



where we assumed D = 2 in evaluating the limits. The diagrams in Figs. 0(c), (d) will require 
the evaluation of the following two limits of l2- 



/ 2 (fc,o,o) 



7-D 



(4tt 

3c 

4n 5 



2 



dt- 



f(l-f) 



[t(l - t)P + ft 2 ] 



7-D 
2 



(?) 4 + o 



oo 



(21) 



J 2 (0,/,0) 



ir 



7-D 



2 (4^)— A) 



dt- 



tf 



3c 
8J> r > 



§(?) a + 0((? 



(22) 



where the limits £ — > oo and p — * are taken with fixed p and ft, respectively, and in D = 2. 



T(± 

We also define c = 



1 

(47r) § 8tt' 

We can now evaluate the series (|Tl~D term-by-term. From the first order term, in the 
limit q — > 0, we get 



Wi(?->0) 



;i + An(o)r 



ft 2 



(23) 



The contributions of order 0(W 2 ) are shown diagrammatically in Fig. |2|. Referring to this 
figure, we label the corresponding terms generated in the expansion accordingly. To illustrate 
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our procedure, we will explicitly calculate the diagram shown in Fig. |2|(e) arising from the 

. is 

^ /dM^^^/i^-^^^y^-^F^on-g-^)), (24) 



final term in ([IT|) . This term is 
A 2 



where 



V(k)V(q- k)V(k')V(-q- k')) = 2 (V(k)V(k')) (V(q- k)V(-q- k')) (25) 
are the contractions which contribute for q^ 0. Using the definition 

W(q)S{r) = (v(q)V(-q + r)) , (26) 



and integrating over k', (gj) becomes 

2A 2 



W 2(e) (q) 



{l + XU(q)Y 
2A 2 



{i + An(o)} 

We now replace W(k) in (P7|) to first order in W to get 



dkh(k, q)h{-k, -q)W{k)W{q- k) 
dkI 2 {k,0)W 2 {k), q^O. 



(27) 



W 2{e) (q^0) 



2X 2 W 2 



dk- 



l 2 (k,0) 



{1 + XH(0)} J {l + \U(k)} 



^ VA 2 ? 



2 A 

Q, 2 / xdx' 

JO 





+ (A)] 




+°(*): 



(28) 



where the last line follows from substituting (19) and (p0| ) into the previous line and making 

k 

the change of variable x = — ; A is the usual ultraviolet cutoff imposed by the lattice. 
Expanding the denominator in (pg|) and integrating over x now yields the result 



W 2(e) (q^0) 



W 



2tt 5 VA 2 c 2 



Q 2 



A\ 2 32 /A\ „ / /A 



(29) 



It is important to note that ( p9|) does not vanish as f2 — > 0. 

The calculation of the remaining terms now follows in a similar way. The diagrams 
arising from the second term on the RHS of (JO]) are those shown in Figs. ^|(a),(b). These 
give 



~ 4 / W 

^«(^0) = -(^ 



si 2 



(30) 
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and 



4tt 3 U 2 c 



W 



16 M 



7T 



+ In 



The third term on the RHS of ([TT| ) gives rise to the diagrams |2|(c),(d). These give 



(31) 



1 



w 



47T 3 VA 2 c 2 



+ 



16 /A 



7T Vfi 



+ In 



(32) 



and 



7T° 



1^ 

AV 



fi 2 



(33) 



Note that the two highest order terms in (^) cancel exactly with those in (|3~TD . Summing 
the contributions (^3|) and (^)-([33|), we then get the result quoted in (^). As mentioned, 
the second order term that remains constant when Q — > comes entirely from the diagram 
0(e). 
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7 Figures 




Figure 1: Diagrammatic representation of Eq. (0). The heavy dashed line represents the 
self-consistently screened random potential, while the thin dashed line is the bare random 
potential. 




(e) 



Figure 2: Diagrams corresponding to the second order terms in the expansion fllT]) . 
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Figure 3: Finite size scaling attempt of the ground state energy e with z = 0.9 demonstrating 
the lack of a transition in our numerical calculations. The disorder averaging was done over 
500 configurations for L = 6, 1200 for L = 8, 1000 for L = 10 and 1000 for L = 12. 
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0.93 




0.85 




Figure 4: Finite size scaling attempt of the ground state participation ratio with Df — 
0.5. Again, the inability to cross these curves at a common point indicates the lack of the 
transition. Disorder averaging was done over the same configurations as in Fig. || 
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